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ABSTRACT 

Using 1-loop renormalisation group equations, we analyze the effect of 
randomness on multi-critical unitary minimal conformal models. We study 
the case of two randomly coupled Mp models and found that they flow in 
two decoupled Mp_i models, in the infra-red limit. This result is then extend 
to the case with M randomly coupled Mp models, which will flow toward M 
decoupled Mp_i. 
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1 Introduction 



Recently, many theoretical models have been proposed in studying the effect of random- 
ness in two-dimensional systems (see for example 0). In a conformal invariant pure 
model, quenched impurities introduce a perturbation term which generally breaks confor- 
mal invariance. The effect of a weak disorder can be anticipated by the Harris criterion 
0. The relevance of disorder in the critical regime of a system can be evaluated by power 
counting. The naive dimensional analysis of the impurities induced term perturbing the 
(conformal invariant) Hamiltonian of the pure system give us the importance of disorder 
near the transition point. One of the first results in the case of uncorrelated quenched 
disorder was obtained in ^, |^, |^ in the context of the Ising model with random bonds 
and in for the Baxter model. These cases correspond to a marginal perturbation of 
the original Hamiltonian. The problem becomes more complicated when quenched im- 
purities induce a relevant perturbation term. Such an example is found in the g-state 
Potts model with random bonds 0, |^ flOl- Using the replica method, it was argued 



that the system reach a non trivial fixed point giving new critical g-dependent exponents. 



More recently, in [Q, 11, 12, 13| it was shown that in some cases impurities can drift a 
first order phase transition to a continuous one, which is in many cases of Ising-like type. 
It was pointed out that, for the g-state Potts model, the g-dependence of exponents for 
different kinds of disorder is not obvious. The absence of non-perturbative results for 
uncorrelated disorder forbid us to give a definitive answer. 

On the other hand, there exist powerful method in studying conformal field theories 
perturbed by some particular operators [Q. Some exact results have been obtained 
by perturbing unitary minimal models. These models corresponds to multi-critical two- 
dimensional statistical models [jl5[, [|1^. It is interesting to study these models when 
the perturbation corresponds to an addition of disorder. This is what we will do in this 
paper. The model that we study consists in taking two copies of such models and add 
a random energy density depending coupling between them. Our result is that under 
particular conditions, our systems will have a large scale behavior corresponding to that 



found in [|1J], namely they flow, in the infrared limit, in a new minimal model. In section 



2 



(2) we show how disorder is implemented in our initial pure systems. In the context of 
the replica method we show which perturbing term we have to consider in the effective 
Hamiltonian of our replicated system. Then, in section (3), we specialize in the case 
of two randomly coupled minimal models. We study the 1-loop Renormalisation Group 
(RG) behavior of the system for an arbitrary number of replica and the quenched case is 
developed in section (4). In section (5) we present a generalization to the case with M 
randomly coupled Mp models with M > 2 arbitrary. Section (6) is devoted to discussions 
and conclusions. 



2 Minimal models perturbed by randomness 

We are interested in conformal invariant multi-critical models perturbed by randomness. 
Our pure system will then consist in unitary minimal models Mp of central charge c = 



^ ~ p(p+i) I — I' I— !■ "^^^ conformal dimension of any operator in the conformal grid of 



the Mp model is A„/ „ = ^^^(p+y (physical dimension is 2/\n',n)- For the unitary 
minimal models that we consider here, the energy operator correspond to e = 02,1 and 
has conformal dimension Ao i = 7 — 4 with e = 

In this letter, we will consider a model consisting in 2 minimal models Mp randomly 
coupled in the following way: if i^o.i and Hq 2 are respectively the Hamiltonians of the 
two pure models, our total Hamiltonian is: 

H = HoA + Ho,2 + Jd^x q{x) ei{x)e2{x) (2.1) 

where q{x) is a random coupling term between the two models. The case of a non-random 
q{x) was already shown to be integrable and leads to a massive theory ||T^. In the context 
of the replica method, we take the average of the n*^ power of the partition function: 

Z^= [ Y[dP[q{x)]Z- (2.2) 

where dP[q{x)] is a symmetric normalized probability distribution for q{x). (|2.2| ) gives 
naively the effective Hamiltonian: 

n n „ 

H = E(^o,i + ^0,2) - ^ E d'x ie\ele{ei) (x) + ■ ■ • (2.3) 

■ 1"^ 



Here we just have write the second cumulant of the distribution P[q{x)]. If p is large 
enough, the dimension of the energy operator will be close to | and naively, higher 
cumulant terms will be irrelevant operators. The quenched case will be obtained in the 
limit n 0. As a first stage we will study the model ( p. 3D for generic n specializing then 
in the quenched case. 

We still have to make an important remark, before going to the detailed calculations. 
We can see that in both cases, the operator algebra (OA) of the energy operator contains 
the operator 03 i and his descendants : 

ee ^ [/] + [03,i] 

The conformal dimension of (p^^i is given by A^^i = 1 — y and so it is a relevant operator. 
In fact, in the interaction terms displayed above, terms with same replica indices or 
coming from higher cumulants of the probability distribution will produce, apart from 
trivial or irrelevant contributions, a term of the form: 



i=l 

Here we have denoted by $ the 03 i operator. The problem of one Mp model perturbed 
by the 03 i operator has been studied extensively in [|I^. There is a non trivial infra-red 
(IR) fixed point, and it has been shown that the system flows to the Mp_i model. This 
fixed point will also be present in the (RG) behavior of our particular model and will be 
of particular interest in the quenched case. 



3 Generic case 

We first consider the case where we have n replicas and thus our model consists in 2n 
minimal models coupled together. The idea is to study the (RG) behavior of ( |2.3|) for 
n > 2 and to identify the different fixed points that can appear. As explained in the 
previous section, the model that we consider in the following is described by the more 
general Hamiltonian : 

H = T.Ki + K2 + d^xY.[^\{x) + %{x)) +gj£xY. {x) (3.1) 

i=l i=l i^j 
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where Hq ^ and Hq 2 are the Hamiltonians of the unperturbed systems, each of them 
corresponding to a minimal model Mp. Note that we have replaced cr by —g, thus the 
physical case for a random model corresponds to g < 0. The 1-loop RG equations can 
be easily obtained from the operator algebra of the perturbing fields : 

n n n 

i^j k^l ij^j 

n 

+4{n - l)Cf,\x - E + ^M) + ■ ■ ■ 

f:($u^)+*2(^)) EH(i/)+<^>2(i/)) - c'i^i^-i/r'^^E(<^>ui/)+<^>2(i/))+--- 

l=\ i=\ i=\ 

n n n 

1=1 i^j i+i 

(3.2) 

where we have omitted the descendent terms. and C%,^ are the structure constants 
of the model Mp. They are symmetric under permutation of the three indices and their 



values can be obtained from 16 



= ^ + 0(e) ; C%^ = ^^ 0(e) 
Using the formula of the 1-loop RG equation llj], [|14| : 



5(i = (2 - dim{gi))gi - TtK'j^gjgk 

(where gi = g, A) and the values of the generalized structure constants Kjg^ obtained from 
eq. ( |3.2| ), we get the following system of equations : 

g = 2eg-{n- 2)g^ ~ Xg + ■ ■ ■ 

A = |A-^-^(n-l)/ + ... (3.3) 
after the trivial redefinitions g g/An; \ X/Air^/S. The first step in the study of the 



RG flow is to find all the fixed points of eq. ( |3.3| ). This is easily done by solving eq. (|3.3| ) 
with the conditions g = X = 0. For > 2 we get 4 solutions : 

A = ; g = (3.4) 
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A = 4e ; ^ = (3.5) 

2e^f2 



2{n - 2)2 + 9{n - 1) 



A = 2e 1 + 



^= , ; A = 2e(l- , ' ' ) (3.6) 

/2(n - 2)2 + 9(n - 1) J2(n - 2)2 + 9(n - 1) 



v^(r2 


-2) 


v/2(^ - 2)2 


+ 9(72-1) 


^/2{n 


-2) 



Note that (|3.4|) and ( |3.5|) are the two fixed points present in the work of Zamolodchikov 
who found that a pure Mp model, perturbed by a 03^i operator fiows to the point 
( ^■51 ) which correspond to the Mp_i model. The next step is to study the stability of 



each of these fixed points. This is done by linearizing (O) around the solutions given 
above g = g* + 6g ; \ = \* + 6\ and getting a linear system : 

The eigenvalues of the matrix A for each of the cases (|3.4|) to (p.6|) give us information 
about the stability of these points. It is easy to see that for ( |3.4| ) both eigenvalues are 
positives indicating that this fixed point is unstable in all directions, while for ( p.5|) we 
have a stable fixed point in both directions. For ( |3.6| ) and ( |3.6| ) we obtain one real 
negative and one real positive eigenvalue, that is, these points are stable in one direction 
and unstable in the other. Thus, they can be reached only if we fine tune the values of 
g and A to keep our system in the stable line of these points. Studying in detail the fiow 
diagram of ( |3.3| ) we can see that the initial conditions A = ; g ^ will fiow far from 
our fixed points toward either a massive theory or another fixed point which can not be 
seen at this order in perturbation theory. 



4 Quenched system 

We now turn to the case n = which correspond to the quenched system of two randomly 
coupled minimal models. By just putting n = in (|3.3| ) we get our new RG equations: 



g = 2e(7 + 2/-A(7 + --- 
; A2 3 2 

A = yA-- + -,2 + ... (4.1) 



Solutions (O) and (3^) are still valid with the same kind of stability but now there is no 
more fixed point solutions with g (points ( |3.6| ) and ( PT^ ) became complex.) So, ( |3.4| ) 
and ( p.5|) are the only fixed points at this order in perturbation theory. Assuming that 
higher loop corrections to ( [4.1|) will no change the qualitative behavior of the fiow near 
our two fixed points, we can see that a system with initial conditions Aq = ; Qq < 
will fiow toward the point ( p.5|) . This is supported by the numerical calculation of the 
RG fiow (|4.1| ) in figure 1 for different values of go. 
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Fig 1 : Flow of the coupling constants for 2 M_p randomly coupledmodels. 

Initial conditions Aq = ; go < are precisely what we expect for the case of a quenched 
system, since the term in the bare Hamiltonian proportional to the $ operator will have 



a factor n, coming from the contraction of the other pairs of energy operators in ( |2.3| ). 
Then, in the limit n = 0, Aq cancels and go = —o < 0. So, adding a small random 
coupling trough the energy densities of two Mp minimal model will drive our system at 
large scales to two decoupled Mp_i models. 



5 Generalization to M coupled models 



In this section we will consider the generalization to M coupled models. The hamiltonian 
is thus replaced by 

M n M n M n 

H=Y.T.H'o,a + \J d'xY.T.{^a{^))+Pj E E (^) (5-1) 

a=l i=l a=l 1=1 a^b,c^d ij^j 

In the following computations, it will be more convenient to express the last term like 

„ M n 

pj £x E i^yb^^d) (x) = 

a^b,c^d i^j 

M M n M n 

9jd'xY.T.T.{^lAeiei){x) + a J d'x E (^44^^) (^) (5-2) 

ay^b c^d i^j <a,b,c,d> i=l 

Here < a,b,c,d > means all the summation over a, b, c, d which take different values 
two by two. Under the renormalisation group transformations, each of these terms is 
going to behave differently and thus there will be two different equations for g and a. 
In fact this is only true for M > 4. For M < 4, the a term is absent. From the set of 
equations with g^ a and A, we will still be able to recover the case M = 2 and M = 3 by 
simply suppressing a. The case M = 2 was already considered in the previous section 
and the case M = 3 will be mentioned at the end of this section. The renormalisation 
group equations for the parameters g, a and A are trivial to compute by generalizing 
computations of previous sections. We obtain 

■g = 2eg-^g^{M{M-l){n-2)+4:{M-2f)-ga{M-2){M-3)-Xg + --- 
a = 2ea-^a2(M-4)(M-5)-^/(n-l)M(M-l)-AcT + --- (5.3) 
A = jX-^X^ -^g\n-l)M{M -If -^a\M -1){M -2){M -3) + ■■■ 
Gomg directly to the quenched case, we found 

g = 2eg- g'^{2{M -2^ - M{M ga{M -2){M -3)- Xg + --- 

a = 2ea-^a\M-A){M-5) + ^g^M{M-l)-Xa + --- (5.4) 

A = ^A-iA2 + ^/M(M-l)2-ia2(M-l)(M-2)(M-3) + --- 
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The next step consists in computing the fixed points associated with these equations. 
Using an e expansion, we got the following points: first, because g — g{...), one set of 
solutions is given by 5^ = and, after some additional computations 



(7 = ; 

(7 = ; 

a = ex ; 

(7 = —ex ; 

where we have defined the following quantities 
4 



A = 
A = 4e 

X = 2e{l-y) 
A = 2e(l+y) 



(5.5) 
(5.6) 
(5.7) 
(5.8) 



X — 



y = 



3x 



(M - 1)(M - 2)(M - 3) + 3(M - 4)^(M - 5)^ 
(M-4)(M-5) 



In addition, there is a second set of solutions with g (with some extra conditions on 
M, see bellow.) After some more tedious computations, we get 



g = -eXj^ 
g = -eX^ 

9 = e^M 
9 = e^M 



(7 — —eZ^X^ 



a — —eZj^Xf^ 



(7 — eZ'ljXj^ 



(7 — eZj^Xj^ 



(5.9) 
(5.10) 
(5.11) 
(5.12) 



and A = 2e - g{2{M - 2f - M{M - 1)) - a{M - 2)(M - 3). We also used the following 
definitions 



— 



-(2(M - 2)2 - M(M - 1)) ± 



M 



X 



M 



2((M - 2)(M - 3) - |(M - 4)(M - 5)) 

^(^(^^(M - 2)(M - 3) + (2(M - 2)2 - M(M - 1)))^ 

+ \{Z^f{M - 1)(M - 2)(M - 3) - ^M{M - 1)^)-^ 



(5.13) 



/m = - 17M3 + 65M2 - 64M + 16 



These last four solutions do not exist for every M. We have some additional constraints: 
/a/ is positive, and the fixed points real, only for M > 13; X^J^ is real only for M < 66. 
Thus the four solutions ( |5 .91 - 15. 12]) are real for 13 < M < 66 while for M > 67 only 
solutions ( |5.10| , pT2D do exist. The next step is to study the stability of these solutions. 
Let first note that the physical initial conditions are 



A = 



a < 



(5.14) 



Because g = g{---), we have the condition g < 0. Then, the last two solutions ( |5.11|) and 
( p.l2| ) can be immediately discarded. Thus only 6 solutions can be relevant. The first 
one {g = a = X = 0) is obviously unstable and can also be discarded. The second one 
((? = cr = 0, A = 4e) turns out to be stable {6g = —2e6g, da = —2e6a, 6X = — |e5A ) and 
this independently of M. This point is just the stable fixed point ( |3.5| ) of the previous 
section. The next two solutions {g = 0,a = ±ex, A = 2e(l =F y)) are also unstable. This 
can be seen by noticing that 



'5a 
6X 



eA. 



6a 
6X 



(5.15) 



with 

and then 

det A± 



2 



3(M-4)(M-5) 2 
(M- l)(M-2)(M-3) -2(M-4)(M-5) 



(5.16) 



--x2(6(M - 4)2(M - 5)2 + 2(M - 1)(M - 2)(M - 3)) < 



(5.17) 



for all M > 4, which implies that there is two eigenvalues with opposite sign. Then we 
remain with the last two solutions, ( |5.9| ) and ( |5.1CI| ) for the previously mentioned values 
of M. Again these two solutions will be unstable. Here, it is not possible to give an 
analytical proof for all M. So we computed numerically the eigenvalues of the matrix 
A± defined by 



6g 
5a 



eA 



± 



6g 
6a 
^6Xj 



(5.18) 
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We found that for each 13 < M < 1000, at least one of the eigenvalue of A_ is positive, 
the same being true for for each 13 < M < 66. 

Finally the last step is to check if under the renormalisation group, we are going 
to reach the only remaining fixed point, ( |5.6|) . Again this was done by constructing 
numerically the flow diagram. We construct it for every value of M between 4 and 1000 
and for each of these flows, we found that we reach the stable fixed point ( ^.61) . 

Before going to the discussion, let's also mention the case M = 3. Renormalisation 
group equations for this case are obtained from (|^) by suppressing the a terms. We 



found 2 solutions to these equations: the trivially unstable solution {g = X = 0) and 
the stable solution {g = 0, X = 4e). Again, we remain with only one solution which 
corresponds to 3 decoupled Mp_i models. 



6 Discussion 

In this paper, we have considered a very particular way of adding randomness for more 
general systems than the well studied Ising or Potts models with random bonds. We con- 
sidered the case of two minimal Mp models randomly coupled and our 1-loop calculation 
shows that (weak) randomness can be easily chosen such that criticality is maintained 
and our system will behave at large distances as two decoupled Mp_i models. Then we 
generalized the study to M randomly coupled Mp. Again we found that at the 1-loop 
order, these models will behave at large distance like M decoupled Mp_i models. In fact, 
the operator algebra of the perturbation induced by randomness produce in the effective 
action a supplementary term which drift our system to a unitary I.R. fixed point. These 
results seem to tell us that the critical behavior of some two dimensional systems in the 
presence of randomness depend crucially on the particular model and the kind of ran- 
domness we are considering, in contrast to what seems to happens in the cases studied in 



0, |Tl|, |T^. We expect that higher loop corrections in our renormalization group equations 
shouldn't modify the qualitative behavior of the flow. However, an analytic solution for 
the coupling flow should make more concrete these conclusion. 
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